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New forms of non-relativistic and relativistic hydrodynamic 
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After a brief account of the derivation of the first-order relativistic hydrodynamic equa- 
tion as a construction of the invariant manifold of relativistic Boltzmann equation, we give 
a sketch of derivation of the second-order hydrodynamic equation (extended thermodynam- 
ics) both in the nonrelativistic and relativistic cases. We show that the resultant equation 
suggests a novel ansatz for the functional form to be used in Grad moment method, which 
turns out to give the same expressions for the transport coefficients as those given in the 
Chapman-Enskog theory as well as the novel expressions for the relaxation times and lenghts 
allowing natural physical interpretaion. 



(N 



(N 



O 
=5 



§1. Introduction 



^ . The relativistic hydrodynamic equation is widely and successfully used in nuclear 

and astrophysics: The dynamical evolution of the hot and/or dense QCD matter 
produced in the Relativistic Heavy Ion Collider (RHIC) experiments can be well 
described by relativistic hydrodynamic simulations. 1 )' 2 ) The suggestion that the 
created matter may have only a tiny viscosity prompted an interest in the origin of 
■ the viscosity in the created matter to be described using the relativistic quantum field 

theory and also the dissipative hydrodynamic equations. The relativistic dissipative 
hydrodynamic equation has been also applied to various high-energy astrophysical 
phenomena, 3 ^ e.g., the accelerated expansion of the universe caused by the bulk 
viscosity of dark matter and/or dark energy. 4 ) We must, however, ask if the theory 
^ of relativistic hydrodynamics for viscous fluids is fully established? The answer is, 

unfortunately, no. 

The fundamental problems on the relativistic dissipative hydrodynamic equa- 
tion may be summarized as follows: (a) ambiguities in the definition of the flow 
velocity 5 )' 6 ) , (b) unphysical instabilities of the equilibrium state, 7 ) and (c) lack of 
causality. 8 ) _1 °) 

Here, we note that the non-equilibrium process may evolve through some stages 
of hierarchical dynamics: 11 ) In the beginning of the time-evolution of an isolated 
prepared state, the whole dynamical evolution of the system will be governed by 
Hamiltonian dynamics that is time-reversal invariant. As time goes on, the dy- 
namics is relaxed into the kinetic regime where the time-evolution of the system 
is well described by kinetic equations which describes a coarse-grained slower dy- 
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namics: The Boltzmann equation for the one-body distribution function is one of 
them; 11 ) Usually the original time-reversal invariance is lost in the description by 
the kinetic equation through the coarse-graining. As the system is further relaxed, 
the time evolution will be described in terms of the hydrodynamic quantities, i.e., 
the flow velocity, particle-number density, and local temperature. In this sense, the 
hydrodynamics is the far-infrared asymptotic dynamics of the kinetic equation. 

In this report, we take the relativistic Boltzmann equation (RBE) 8 ) as a typical 
kinetic equation and explore the basic problems with the relativistic hydrodynam- 
ics. For obtaining the proper relativistic hydrodynamic equation, it is a legitimate 
and natural way to start with RBE which is Lorentz invariant and does not have 
stability nor causality problems. 8 ) For analyzing the problems (a) and (b) first, we 
derive hydrodynamic equations 12 )' 13 ) for a viscous fluid from RBE on the basis of 
a mechanical reduction theory called the renormalization-group(RG) method 14 ) -16 -* 
and a natural ansatz on the origin of dissipation. We then proceed to the causality 
problem by deriving the so called extended thermodynamics: 17 )~ 19 ) We shall de- 
rive the mesoscopic dynamics of the RBE by constructing the invariant /attractive 
manifold incorporating the first excited as well as the zero modes of the linearized 
collision operator. We obtain the same expressions for the transport coefficients as 
given by the Chapman-Enskog method 20 ) and also the relaxation times in terms of 
relaxation functions which have natural physical interpretation. We also show that 
the resultant equation suggests a novel ansatz for the functional form to be used in 
the moment method. 21 ) 



§2. Relativistic Boltzmann equation 

The relativistic Boltzmann equation reads 8 ) 

P ft d ll f p (x) = C\f] p (x), (2-1) 

where f p (x) denotes the one-particle distribution function defined in the phase space 
(x , p); p^ denotes the four- momentum of the on-shell particle, i.e., p^p^ = p 2 = m 2 
and p° > 0. The r.h.s. is the collision integral 

C[f] P (x) = 2\ J2 ~o E~o 5^-0 w (P.Pl|P2>P3)(/p 2 ( a; )/p3( a; ) ~ fp( X )f P i( X )}^ 2 ) 

pi 1 P2 ^ p 3 P3 

where uj(p , p\\pi , ps) denotes the transition probability due to the microscopic two- 
particle interaction. Owing to the symmetry property u(p,Pi\p2,P3) = ° J {P2,P3\p,Pi) 
= uj(pi,p\ps,p2) = w(j)3,p2\pi,p) and the energy- momentum conservation, the colli- 
sion operator satisfies the following identity for an arbitrary vector ip p (x), 

P P PI P2 P'A 

Xw(p, Pl\p2 , Pz) (Vp(x) + <Ppi(x) ~ PP2( X ) ~ Ppsix)) 

x(f P 2(x)fp 3 (x) - f P (x) f Pl (x)) , (2-3) 
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from which one can show that the function <p p (x) = a{x) + p^b^(x) is a collision 
invariant satisfying the equation ^ p ^ <Pp{x) C[f] p (x) = 0, where a(x) and b^{x) 
being arbitrary functions of x. This form is, in fact, the most general form of a 
collision invar iant. 8 ) 



§3. Reduction to hydrodynamic equation 



Before developing our analysis based on the RG method, we briefly summarize 
ad hoc aspects in the standard methods such as the Chapman-Enskog expansion and 
Maxwell-Grad moment methods. 8 ) 

In Chapman-Enskog method, 20 ) the zeroth-order solution is given by the lo- 
cal equilibrium distribution function, i.e., the Jiittner function 22 ) fp°\x) = fp q (x). 
Then, one imposes a condition of fit or matching condition, which usually consists 
of drastic assumptions that internal energy and the particle-number density in the 
nonequilibrium state are the same as those in the local equilibrium state. 8 ) Indeed 
these constraints are equivalent to a physical assumption that there is no internal 
energy nor particle-number density of dissipative origin. Such an assumption, how- 
ever, does not have any solid physical foundation, since the distribution function in 
the nonequilibrium state should be quite different from that in the local equilibrium 
state. Such matching conditions are also imposed in the moment method. 8 ) 

In the RG method that we adopt, 15 )' 16 ) one needs no such matching conditions 
for derivation, and their correct forms are obtained as a property of the derived 
equation once the frame is specified by the macroscopic-frame vector. 12 )' 13 ) We 
have shown 12 )' 13 ) that the conditions of fit in the energy frame are compatible with 
the underlying Boltzmann equation, but those in the particle frame adopted in the 
literature are not and thus will never be satisfied in any physical system. 

We solve RBE (2-1) in the hydrodynamic regime, and thereby derive the hy- 
drodynamic equations governing the hydrodynamic variables which are introduced 
to parametrize the distribution function. To make a coarse graining with Lorentz 
covariance, 12 )' 13 ) we introduce a timelike Lorentz vector a M , with a > 0. Thus, 
specifies the covariant and macroscopic coordinate system where the velocity field 
of the hydrodynamic flow is defined. We call a M the macroscopic-frame vector. Al- 
though a M could depend on the momentum p as well as the space-time coordinate 
x, i.e., = a,p(x), the possible momentum dependence may not be legitimate as a 
macroscopic frame vector*). Once is given, RBE(2-1) is rewritten in terms of the 
new coordinates (r , a^) as 

, a) = ) C[f\{r , a) - e p ■ V/ P (r , a), (3-1) 

or p- a p (T , a) F p- a p (T , a) 



*-* As is shown in 12), 13), the Eckart frame can be realized only when the macroscopic frame 
vector has a p dependence. This implies that the macroscopic space-time in the Eckart frame is 
defined for respective particle state with a definite energy-momentum, which may have a difficulty in 
a physical interpretation. 23 - 1 In fact, it can be shown 23 - 1 that the p-independent timelike macroscopic 
frame vector with the Lorentz covariance uniquely leads to the hydrodynamics in the Landau-Lifshitz 
frame. 
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where d/dr = (I / a 2 p {x)) a v p (x) d v and V = A p lv {x)d v = d/da^, with A p w {x) = 
g^ v — a p (x) a p (x)/a p (x). In Eq.(3-1), the small parameter e represents the space- 
nonuniformity, which may be identified with the ratio of the average particle distance 
over the mean free path, i.e., the Knudsen number; e will be set back to unity in the 
final stage. This seemingly trivial rewrite of the equation expresses the assumption 
that only the spatial inhomogeneity is the origin of the dissipation. It is noteworthy 
that the RG method applied to the nonrelativistic Boltzmann equation with the 
corresponding assumption successfully leads to the Navier-Stokes equation; 24 ) the 
present approach 12 )' 13 ) is simply a covariantization of the nonrelativistic case. 

In the RG method, 12 )' 13 )' 15 )' 16 ) we first try to obtain the perturbative solution f p 
to Eq. (3-1) around the arbitrary initial time t = tq with the initial value / p (to, c); 
f p {r = To , a ; To) = / p (to , <r). Here we have made explicit that the solution depends 
on the initial time To; we suppose that the initial value is on an exact solution. 
The initial value as well as the solution is expanded with respect to e as follows: 
fp(r, a; tq) = fp°\r, a; To) + ef^\r, a; tq) + • • • , and an obviously similar form for 
fp(r ,a). 

It is easy to see that the zeroth-order equation reads 

(r, a- to) = ) rC[/< >]p(T, a- r ). (3-2) 

or p p- a p (a ; r ) 

Since we are interested in the slow motion that would be realized asymptotically 
as r — > oo, we should take the following stationary solution or the fixed point, 
&pfp°\ T i a 'i T o) = 0) which is realized when C[/^°- ) ] p (t, a; To) = 0, for arbitrary 
a. Thus, we find that the zero-th order solution fp°\r , a ; To) is a local equilibrium 
distribution function 



/W(r, a; r ) = ^3 exp 



l4p_ \ r ) -p^u^a; T ) 
T(a; to) 



= / p 6 >;t ), (3-3) 



with u M (<7 ; To) u^(a ; To) = 1. Note that the would-be integration constants T{a ; To), 
n(a ; To), and u M (cr ; to) are independent of r but may depend on To as well as a. 
Given the zero-th order solution, the first-order equation reads 

^ 1) W = ^^ ? /W(t) + F p , (3-4) 

where A pq denotes the matrix element of the linearized collision operator A given by 

1 (9 

{A) pq = A pq = — — — TrrC[f] p . (3-5) 
P ' a p u h /=/oq 

and 



F P = -— -p-Vf?. (3-6) 
P u p 
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It is essential for obtaining a slow motion to clarify the spectral properties of A. For 
this purpose, it is found convenient to convert A to another linear operator, 

L = (/ eq ) -l A/ eq j (3.7) 

with the diagonal matrix f P q = f p q 5 pq . To characterize the properties of the linear 
operator L, let us define an inner product between arbitrary nonzero vectors ip and 
ip by 

(ip, = E \(p- a p)fp q( Pp^P- ( 3 - 8 ) 

p P 

Then it can be shown 12 )' 13 ) that the linearized collision operator has a remarkable 
properties that it is semi-negative definite and has five zero modes given by 

f p M for a = a, 

^ 1 f A (3 " 9) 

I lxm for a = 4. 

The functional subspace spanned by the five zero modes is called Po space and 
the projection operator to it is denoted by Po; we define Q$ by Qq = 1 — Pq. 
In the following, we also use the modified projection operators defined by Po = 

/eqp Q (/eq)-l and Q Q = f«lQ (f«l)-\ 

We proceed to solve the perturbed equations up to the second order. Then 
summing up the perturbative solutions, we have an approximate solution around 
r ~ t : f p (r ,a;r ) = # 0) (t , a ; r ) + e ff\r , a ; r ) + e 2 ff\r , a ; r ) + 0(e 3 ). 
It should be noted that this solution contains the secular terms that apparently 
invalidates the perturbative expansion for r away from the initial time To- We can, 
however, utilize the secular terms to obtain an asymptotic solution valid in a global 
domain. 15 )' 16 ^ Indeed we have a family of curves f p (r, a; To) parameterized with 
To: They are all on the exact solution f p (a ; r) at r = To up to 0(e 3 ), although only 
valid for r near To locally. Then, the envelope curve of the family of curves, which 
is in contact with each local solution at r = to, will give a global solution in our 
asymptotic situation, which is shown to be the case. 15 )' 16 ) According to the classical 
theory of envelopes, the envelope that is in contact with any curve in the family at 
t = tq is obtained 15 ) by 



d ff 



= 0. (3-10) 

T =T 

The derivative w.r.t. to hits the hydrodynamic variables, and hence we have the 
evolution equation of them that is identified with the hydrodynamic equation. 12 )' 24 ) 
We also note that the invariant manifold which corresponds to the hydrodynamics in 
the functional space of the distribution function is explicitly obtained in the present 
method. 12 )' 13 ) 

Putting back e = 1, Eq.(3T0) is nicely reduced to the following form, 

( P -a p )^- + ep-V [f?- [A^Q F] p }+O(e 3 ) = 0. (3-11) 

p L J 
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If one uses the identity (p ■ a p ) ^ + p ■ V = Eq.(3-ll) is found to have the 

following form 

<VC = 0, <VVf st = 0. (3-12) 

with 

T M-EV^{/; q -[^Qof] p }, (3-i3) 

v p 

Kt^rn-i ^P^l{f7- [^QoFl}- (3-14) 
v p 

3.1. Example: Landau- Lifshitz frame 

In this subsection, we present the hydrodynamic equation derived by the RG 
method. 

If we take the most natural choice for the macroscopic vector as a^(x) = u p (x), 
the resultant energy-momentum tensor and particle current turn out to be 12 )' 13 -' 

Tj£ = e u>* u v - (p - C V ■ u) A» v + 2 r) A^ vpa V p u a , (3-15) 

N p st =nu p + \±V^, (3-16) 

respectively, with A» v = - u» u u , V = AH" d v , and A^ pa = 1/2 ■ (A^A ua + 
A!"* A vp — 2/3 • A pu A pa ). Here, n, p, and /j denote the particle-number density, 
pressure, and reduced enthalpy per particle, respectively. It is clear that these for- 
mulae completely agree with those proposed by Landau and Lifshitz. 6 ) Indeed, the 
respective dissipative parts 5T pu and 5N 11 in Eqs. (3-15) and (3-16) meet Landau 
and Lifshitz's constraints 

Se = u^ST^u,, = 0, 5n = u^5N p = 0, Q M = A^8T vp u p = 0. (3-17) 

Thus we have derived the dissipative hydrodynamics in the Landau-Lifshitz (energy) 
frame in the RG method. 

Moreover, since our theory is in the level of statistical mechanics, we have the 
microscopic expressions for the transport coefficients appearing in the hydrodynamic 
tensor (3- 15) and current (3-16), as follows: 

( = ~(n,L- 1 n), a = ^(j>,l- 1 j m ), r l = ~{^,L-H iaf ).^\%) 

Here, we have introduced the following microscopic thermal forces (IJ p , J p , tt p u ) = 
(IJ P , Jj}, ■Kp lu )/(p ■ u), where 

n p ^(^- 1 ){p-uf+(( 1 -l)Th- 1 T){p-u)- l -m 2 , (3-19) 

Jp = -{{p-u)-Th) A» v p v and vr^ = A^ pa p p p a . Here, 7 denotes the ratio of the 
heat capacities. 
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It is noteworthy that the transport coefficients can be rewritten in the Green- 
Kubo formula. 11 ^ With the use of the "time-dependent thermal force" defined by 
n p (s) = J2q [ eSL ]pq^-Q an d so on > the relaxation functions are given by the time- 
correlators 

R c (s) = ^{n(Q),n(s)), (3-20) 

and so on for R\(s) and R v (s) with obvious modifications to R^(s). Then the 
transport coefficients given in Eq.(3-18) are rewritten as follows, 

poo poo poo 

C = ds R c (s), A = ds R x (s), r) = ds R v (s). (3-21) 
Jo Jo Jo 

We remark that the instability problem of the local equilibrium state is analyzed 
in 13) and 25). 

§4. Second-order equations and moment method 

In the derivation of the first-order hydrodynamic equations like Landau-Lifshits 
equation, we utilized the zero modes of the linearized collision operator, which form 
the invariant manifold on which hydrodynamics is defined; the would-be constant 
zero modes acquire the time-dependence on the manifold by the RG equation. Our 
formalism can be extended to include excited modes as additional components of 
the invariant/attractive manifold. 26 ) The out-come is nothing but the extended 
thermodynamics or Israel-Stewart type equation, with new microscopic expressions 
of the relaxation times and lengths. Furthermore, we find the proper ansatz for 
the distribution function to be used in the moment method. 26 ) For the shortage of 
space, we here only present the results, leaving the detailed derivation in a separate 
paper. 26 ) We emphasize that our theory gives an explicit construction of the invariant 
manifold corresponding to thirteen moments, which has been long sought for. 18 )' 19 ) 

4.1. A brief review of Grad's thirteen-moment method and Grad-Muller equation: 
non-relativistic case 

In Grad's thirteen-moment method, 8 )' 21 ) the one-particle distribution function 
fv(t, x) is expressed in terms of the equilibrium distribution function f^(t, x) as 

fv(t, x) = f^(t, x) ^1 + <Pv{t, x)^j, where <Pv(t, x) is assumed to have the form 

X) = X) = 7$(i, x)TT ij (t, X) + J l v (t, X) J\t, X). (4-1) 

Here, 7riy(t, x) and Jy(t, x) are defined by 7r l y(t, x) = m(5v l (t, x)Sv : >(t, x) — 
\8 ij \8v(t, x)\ 2 ) and J\,{t, x) = (f \Sv(t, x)\ 2 -^T(t, x))5v i (t, x), with 8v(t, x) = 
v — u(t, x) being the peculiar velocity. 

Then the evolution equation of the thirteen coefficients are determined by the 
equations all of which are derived from the Boltzmann equation with use of the 



linearized collision operator given by L v y, = (/^ q ) 1 gj^C[/]v 



/^ q . Thus the 
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Grad-Muller equation is obtained as a closed system of the equations governing T 
n, u l , tt^ , and J 1 , which includes 

d 



m n 7— u 1 + mnu ■ Vu* 
at 



-V 3 ( P S^ -2 v G TT ji 
X l J + (other terms), 



ir lJ + rl 



(4-2) 
(4-3) 



with p being the pressure, where we have defined the thermodynamic force given by 
X'i j = 1/2 ■ (VV + VV - 2/3 • S ij V • it). 

In the Grad moment theory, the transport coefficients and relaxation times are 
given in terms of the inner product defined by (ip , x) cq = ^2v fi^^vXv- For 
example, 



ri G = - 



1 (^,^) cq (vr M ,^) cq G 



eq 



(4-4) 



10T ( TT mn , L TT mn ) eq ' (7T kl , L7T H ) cq ' 

These formulae are different from those given in the Chapman-Enskog expansion 
method; for instance, r/ CE = — 1/10 T • (ir^ , L~ l 7r lJ : ) cq , which clearly shows that 

1] 



,CE 



4.2. T/ie mesoscopic dynamics from the RG method: nonrelativistic case 

We apply the RG method to derive the extended thermodynamics or mesoscopic 
dynamics 18 )' 19 ) from the nonrelativistic Boltzmann equation, which is equivalent to 
construct the invariant manifold 16 ) in the space of the distribution function. It is 
nice that we can read off the proper ansatz for the deviation function <Pv from the 
distribution function constructed in our method. Leaving the detailed derivation of 
the formulae to a separate paper, 26 ) we here give the results: Our deviation function 
<P V and the relaxation equations read 



<2> 



(4-5) 



and 



V 



9 TT 

dt +V ' V 



fv = Y,Y,fv^[L- l ^] v L vk ^ (4-6) 



« k 



d_ 
di 



+ v ■ V 



fv 



1 



v k 

respectively. The corresponding microscopic expressions of the transport coefficients 
are given by 



V 



TK 



1 



A TK = - 



\J i ,L~ 1 J i ) ( 



(4-8) 



10 T s " ' — " /eq ' " 3 T 2 x " ' " " /cq ' 

which perfectly agree with those by Chapman-Enskog method, r/ CE and A CE . Fur- 
thermore, we have the microscopic representation of the relaxation times given by 
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which are new and different from those in the naive Grad moment method. 

4.3. The mesoscopic dynamics from the RG method: relativistic case 

In the relativistic case, if we express the distribution function by f p {x) = fp q (x) (l+ 
<P p (x)), our RG method gives the following expression of @ p (x), 

P = "f E L vl fa n + J? J, + < TT^), (4-10) 

q 

where Up, Jp, and ftp" are the microscopic thermal forces, This form is different from 
those used by Israel-Stewart 10 ) and Denicol et al. 27 ) within the moment method. 
The relaxation equations derived in our RG method read 

II = — V • u — Tn DII + other terms involving relaxation lengths, 

J M = J- V M — — tj A^ a DJ a + other terms involving relaxation lengths, 
h 2 T 

vr^ = A^p* V p u a - r w A^ ab Dir ab 

+ (k^ A^ V-u + «W A^ ac A^ A ahde V V + 4 2 2 z^ ac 2\ c 6 ^ Wa6 ) vr pa 

+ other terms involving relaxation lengths, (4T1) 

with the volticity oj^ u = \ (V fl u u — V^u^) and D = d^. The energy-momentum 
tensor and particle current read T^ d = e u M u v — (p + C H) ^ v + 2 r\ ■K tML ' and iV^ nd = 
+ A J M , respectively. 

The RG method gives microscopic expressions for the resulting times Tn, tj, 
and tv as follows 

77 ~ (n,L-in)~ f °°dsR c (s)> (4 ' 12) 

and so on for rj and T n with obvious modifications. It is noteworthy that the relax- 
ation times Tn, tj, and T n are represented in terms of the relaxation functions R((s), 
R\(s) and R v (s): These formulae have a clear physical meaning of the relaxation 
times as the correlated times in the relaxation function. We note that it is for the 
first time that such representation and the interpretation of the relaxation times are 
given for the extended thermodynamics, as far as we are aware of. We also note that 
shear viscosity derived in our method do coincide with that in the Chapman-Enskog 
method as shown before. 



§5. Brief summary 



We have reported our attempts to derive first-order and second-order relativistic 
hydrodynamic equations from relativistic Boltzmann equation which has a manifest 
Lorentz invariance and does not show any pathological behavior such as the instabil- 
ity and acausality seen in existing hydrodynamic equations. We have given the novel 
extended thermodynamics both in nonrelativistic and relativistic cases through the 
explicit construction of attractive manifold containing the relaxation process from 
Boltzmann equation. 
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